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We calculate the Casimir force between slabs of finite thickness made of intrinsic and doped silicon 
with different concentration of carriers and compare the results to those obtained for gold slabs. We 
use the Drude and the plasma models to describe the dielectric function for the carriers in doped Si. 
We discuss the possibility of experimentally testing the appropriateness of these models. We also 
investigate the influence of finite thickness on VO2, which has recently been proposed for Casimir 
effect measurements testing the metal-insulator transition. 



I. INTRODUCTION 



The availability of experimental set-ups that allow accurate measurements of surface forces between macroscopic 
objects at submicron separations has recently stimulated a renewed interest in the Casimir effect [Ij and its possible 
applications to micro- and nanotechnology. 

In 1948 H. Casimir calculated the force between two plane-parallel mirrors placed in vacuum at a distance L apart 
from each other and with the area A of the mirrors being much larger than the squared distance (A ^ i^) . In the 
ideal case of perfectly reflecting mirrors the force is given by the following expression 

with a positive value of Fcas corresponding to attraction, and a subsequent negative pressure. 

The Casimir force can be understood as the effect of radiation pressure on the Fabry-Perot cavity formed by the 
two plane and parallel mirrors. The intracavity vacuum energy being either enhanced or suppressed, depending on 
whether the field frequency is resonant or antiresonant, the net Casimir force results from the balance between the 
repulsive and attractive contributions associated respectively with these frequencies. The force is then obtained as an 
integral over the axis of real frequencies, including the contribution of evanescent waves besides that of ordinary waves, 
and transformed into an integral over imaginary frequencies by using physical properties fulfilled by real mirrors used 
in experiments which do not require any further adhoc hypothesis. 

Compared to the ideal situation considered by Casimir, a number of corrections have to be taken into account in 
the calculation of the force in real experiments. A large number of papers have been devoted to the study of these 
effects and we refer the reader to for an extensive bibliography. Here we will only be concerned by the influence of 
the material properties and slab thickness. 

Considerable experimentalprogress has been achieved 0] in the control of the Casimir effect, opening the way to 
applications in nano-science [3l, '4] , particularly in the development of nano- or microelectromechanical devices (NEMS 
or MEMS). NEMS are movable nano mechanical structures inspired from MEMS, with minimal critical dimensions of 
a few tens to a few hundreds nanometers. NEMS approach constitutes a real technological breakthrough to prepare 
the future generations of sensors and actuators. At such small distances between the different elements, the Casimir 
force in these systems may become quite important. It may, on one hand, perturb the systems and produce stiction 
and adhesion [3|, but also be put to good use as an external force allowing to change the systems' resonance frequency 
or introduce bistable behavior [1]. 

In the past 10 years the Casimir effect has been studied extensively for metals. Experiments have been performed 
for Au, Al, Cu in different geometries and experimental set-ups. On the theoretical side, calculations have taken into 
account the finite conductivity and frequency dependent reflection coefficient of the different metals, modeling them 
either by plasma or Drude model or taking into account tabulated optical data[3, H, S 0, Hi Hi- Only recently, a 
Casimir force measurement for Silicon bulk mirrors has been reported on [Tl| and calculations of the Casimir force 
for Si bulks and slabs (1^ il3i] have been performed however without investigating further the observed difference in 
the bulk and slab behavior. The temperature dependence of the force between Si plates has been studied in [14], while 
the influence of skin depth on the Casimir force between metallic surfaces has been observed recently by Capasso and 
collaborators [15.] . 

The reference material in nano- or micro-electromechanical devices is of course Silicon. In this paper, we present 
calculations of the Casimir force between slabs of intrinsic or doped Silicon, Gold and VO2. In particular, we 
concentrate on the influence of the slab thickness on the value of the force. We use the Drude and the plasma model 
for the dielectric function of conductors and a Drude-Lorentz model for Si and the insulating VO2 . We have previously 
shown that for intrinsic Silicon the Casimir force depends strongly on the slab thickness and reduces considerably 
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if the slab separation exceeds the slab thickness, while for Gold force reduction is negligible, except for very thin 
slabs [l^ • If the conductor is modeled by a plasma model (7 = 0) the Casimir force becomes independent of the slab 
thickness. Our present calculations show that for doped Silicon, the Casimir force is diminished only in some distance 
range and for low carrier levels and re-increases at long distances. We explain these effects by analyzing the phase 
factor acquired by the vacuum field due to finite slab thickness and the different behavior of the dielectric functions. 
We also discuss the effect of a Gold coating on the Si slab. Even though, a description of a thin Gold film by a local 
dielectric function is not completely reliable anymore, it gives some hint about the expected behavior. The slightest 
Gold coating on intrinsic Silicon, even of only Inm thickness, re-increases the Casimir force between the two Si slabs. 
The above described decrease in the force due to finite thickness of the Si slab vanished because of the Gold coating. 

We finally analyze the influence of slab thickness for VO2 which has the particularity to undergo a metal-insulator 
transition at 340 K [13, [3. An experiment was proposed recently [l^ to investigate in detail the temperature 
dependence of the Casimir force between gold and VO2 film. The temperature dependence of the Casimir effect is 
still an issue of controversial discussions j2l|, [H, [H, [13, [H, [H, H^, [M, HI, M, HI, [H, [13, [11. Here we make 
account only of the crucial change in the conductivity due to the phase transition, but do not consider the temperature 
corrections to the Casimir effect itself. 



II. FORMULATION OF THE CASIMIR FORCE BETWEEN SLABS OF FINITE THICKNESS 

The Casimir force is usually written as an integral over imaginary frequencies and wavevectors jssj . In order to 
visualize easily the variation of the real Casimir force with respect to the ideal formula ([1]) it is convenient to introduce 
a reduction factor rip = F/Fcas 01 

12{)L^ f°° f°° 

VF^^^J2 / / dc.Ac/(^,fc), (2) 
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Tp denotes the reflection amplitudes of the mirrors at a given polarization p. This notation is a shorthand for Vp [ilu, in) 
where ito is the imaginary frequency and in the imaginary wavevector along the longitudinal direction of the cavity 
while k is the modulus of the transverse wavevector. The second integral is written over imaginary frequencies as 
explained in detail in 

Assuming the plates to have a large optical thickness, the reflection coefficients Vp correspond to the ones of a mere 
vacuum-metal interface 43l 



\J ijj"^ {e {iuj) — 1) + c^K^ — CK 
■\/ Lo"^ (e {iiS) — 1) -I- c^K? — CKE (iuj) , , 

P\\ — , , V*^/ 

\/uj'^ (e (iuj) — 1) + c^K^ + CKE {iuj) 

Pp stands for pp (iLo^in) and e {iuj) is the dielectric constant of the metal evaluated for imaginary frequencies. 

However, for thin mirrors or slabs the reflection coefficients depend on the physical thickness D and evaluate to 
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D 



c 



y/uj^ (e {iij)-l)+C^K\ (4) 



£ (iLu) is the dielectric function of the material evaluated for imaginary frequencies, d is the optical length of the slab. 
The single interface expression is recovered in the limit of a large optical thickness S ^ 1. 

The dielectric response function for real frequencies may be written in terms of real and imaginary parts e' and 
e" obeying usual causality relations which allow one to obtain the dielectric constant at imaginary frequencies e (iuj) 
from the function e" (x) evaluated at real frequencies x [i^ 

ei^u;)-l^- H ^f^dx. (5) 

TT In X^ -f 
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FIG. 1: The difFerent dielectric functions for intrinsic and doped Silicon for varying carrier densities in comparison with Gold. 
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TABLE I: The values of plasma frequency and relaxation rate for various carrier densities, uip — Ne^ / {eom*), 7 = Ne^ p/m* , 
where m* — 0.34me is the effective mass of the holes, and p is the resistivity [l3l ]. 



The tabulated optical data for the complex index of refraction for Silicon and Gold can be found in [37|, while the 
data for VO2 has been given in a dedicated paper For Silicon and VO2 the data covers the whole relevant 

frequency range and, in contrast to most metals, no extrapolation procedure at low frequencies is necessary. We 
will give frequencies either in eV or in rad/s, using the equivalence 1 eV = 1.519 x 10^^ rad/s. Figure [1] shows the 
dielectric permeability of doped Silicon with difFerent carrier concentrations as a function of imaginary frequencies, 
obtained by using ([5]), which we will need for the calculation of the reduction factor of the Casimir force While 
at low-frequencies the dielectric function of intrinsic Silicon approaches a constant value — 11.87, e{uj) of doped Si 
increases with increasing carrier concentration and behaves like that of diluted metal. With increasing frequency, the 
dielectric function of intrinsic Silicon is nearly constant up to about 10^^ rad/s and falls off only for high frequencies 
above a cut-off frequency, ujq k, 6.6 • 10^^ rad/s, towards its asymptotic value £00 = 1.035. £(w) of doped Silicon 
decreases rapidly with increasing frequency and shows the same cut-off frequency as intrinsic Si. 

The dielectric function of intrinsic Silicon can well be approximated by the following Drude-Lorentz function (4^ 

esi[iuj) = £00 H ^ — (6) 

with the numerical values as given above. The dielectric function of doped Si contains an additional part, which is 
modeled with a dielectric function given by a Drude model 

edop(«w) = £si(«cj) -f — - — ■ — -. (7) 
ijj{uo -1-7) 

The different values of the plasma frequency Wp and the relaxation rate 7 for various carrier densities are given in 
Table H 

The description of the dielectric function of the p-doped Silicon by the model ([7]) is valid for the doping levels up 
to 10^° cm~^. For higher doping levels Si becomes degenerated. 

The dielectric function of VO2 is shown in Fig. above and below the critical temperature T^. Below the critical 
temperature we used the model already applied in [l9| which had been proposed first in [l7j 

I- \ 1 I £n(*Oo) — 1 ^ S„ j 

erAiLo) = 1 H 2 \- / 2 (8) 
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TABLE II: The values of the parameters in the models (O and ^ [ij 




FIG. 2: Dielectric function for VO2 above and below the critical temperature in comparison with the one of intrinsic Si, Gold 
and AI2O3. 



with e„(ioo) = 4.26, uJoo — l5eV. For the rest of the parameters see Table HIl Above the critical temperature the 
dielectric permeability is given by 



eft(zw) = 1 + 



uj{uj + 7fi) 



^— 1 



(9) 



where eft(icx)) = 3.95, ojp^n = 3.33 eV, 7,1 = 0.66 eV. 

Clearly, above the critical temperature the material behaves more like a metal with a strongly increasing dielectric 
function at low frequencies. For low temperatures its dielectric function nearly matches the one of intrinsic Si, as VO2 
becomes an insulator and the dielectric function at zero frequency has a constant value of about e„(iO) ^ 10. 



III. NUMERICAL RESULTS FOR DOPED SILICON 



We now evaluate the Casimir force ([2]) using^he tabulated optical data for different slab thickness D. The numerical 
procedures follow the principles described in [5|. 

Figures [3] and [4] show 77 as a function of plate separation for intrinsic and p-doped Silicon for bulk mirrors and for 
Si slabs of lOOnm thickness. Between two bulk mirrors the force reduction factor is a continuously growing function 
of the mirror separation and reaches a constant value r]p{oo) ~ 0.303 in the long distance limit, which means that 
F « -Fcas/3. In contrast, between two Silicon slabs the Casimir force reduction factor grows continuously only for 
separations L ~ D, while it starts to diminish considerably when the slab separation becomes of the order of or 
exceeds the physical slab thickness [l^. This effect disappears for doped Silicon. Only a small decrease in the force 
reduction factor remains at intermediate distances (about 1/im). For larger distances the conduction due to the 
carriers introduced by doping re-increases the value of the reduction factor. 
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FIG. 3: Reduction factor of the Casimir force between two bulk mirrors of intrinsic or p-doped Silicon. 




FIG. 4: Reduction factor of the Casimir force between two lOOnm thin slabs of intrinsic or p-doped Silicon. 



This result could have very interesting consequences in nanotechnology. The observed decrease in the Casimir force 
for intrinsic Si in the distance range L > D can be canceled by injecting carriers. This allows to tune the Casimir force 
within certain restrictions using at the same time the slab thickness and the carrier densities as variable parameters. 

In [llj the first experiment of optical modulation of dispersion forces through the change of carrier density by laser 
pulses was reported. The irradiation of a silicon slab by laser pulses allows to achieve charge carrier concentrations 
of n = (2.0 ± 0.4) X 10^^ cm~"^. Here we present numerical results for the Casimir force obtained by assuming that 
the contribution of the induced charge carriers is given by Drude or plasma terms, oj^ ^/ {uj{uj + 7^)) or oj^ il^'^ ^ 
i = e,p, with LUp^p = 0.368 eV, jp = 0.00329 eV for the positive induced carriers and uip^e — 0.329 eV, je = 0.01185 
eV for the negative ones. 

Fig. [5] gives the numerical results for the reduction factor of Casimir force between two silicon slabs irradiated by 



laser pulses in comparison with p-doped silicon. Here we took the parameters from the paper IjJ. The doping is 
described by the Drude or plasma model. The difference between the models could manifest itself at distances of 
several micrometers for highly doped silicon provided the slab is thin (upper graph) . However the experiment requires 
a mirror thick enough for irradiation not to penetrate into the cavity. The thickness should be greater than the optical 
absorption depth of Si at the wavelength of the laser beam used for the doping. If the thickness of the slab is 4000 
nm as used in flTl the experiment can hardly distinguish between Drude and plasma models. 



IV. NUMERICAL RESULTS FOR VO2 

In [1^ a compound mirror, having a lOOnm layer of VO2 on an AI2O3 bulk substrate, is treated as an effective 
medium with a dielectric permittivity (j8l9p . The authors rely on the fact that the parameters in Table IIIl were retrieved 
from reflectivity and transmissivity spectra exactly for the same two-layered system as in Ref. [17]. 

However, [l7[ contains also the description of the fitting procedure that uses the formulas corresponding to two 
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layered system. The finite thickness of the substrate is neglected in the analysis, but the account is made for the finite 
thickness of the VO2 film. The dielectric function of the substrate (sapphire) is considered as known. It is substituted 
into the formulas, and the unknown parameters of the model describing the VO2 film are found by the curve fitting 
programme. Since it is not quite clear from T?] if the parameters in the Table HIl and the model correspond to the 
effective media or to the VO2 itself, we considered both possibilities. 

First we calculated the reduction factor for the Casimir force between two "bulks" of the effective medium. In 
other words we use the dielectric permittivities ([5]) and ^ with the parameters from Table |TT] and formulas for bulk 
refiection coefficients. The corresponding curves in Fig. [S] are plotted by thin full line below the phase transition, 
T < Tt, and by thick full line above it, T > T*. 

The reduction factor for the Casimir force between two 100 nm VO2 slabs below the phase transition is given by 
thin dotted line. We see that before the phase transition the behavior of VO2 slabs resembles the one of the silicon 
slabs. Above the phase transition (T > Tj, thick dotted line) the material behaves like a dilute metal, and we do not 
observe the decrease of the reduction factor at L > 200nm. However the reduction factor for a slab never reaches the 
unit value it does for pure metals at zero temperature in the long distance limit. 

We also calculate the reduction factor for lOOnm VO2 film on a sapphire substrate. It means that we used the 
formulae for two-layered mirrors. The layer facing the cavity has the dielectric permittivity ([5]) at T < and © at 
T > Tt- For the substrate we used the dielectric permittivity of AI2O3 

with Ai = 1.023, A2 = 1.058264, A3 = 5.280792, and /i = 20.19 eV, /2 = 11.21 eV, h = 0.07 eV [H. 

Below the transition temperature the dielectric function of VO2 is close to the one of AI2O3 (see Fig. [5]). That 
is why the result for the Casimir force between two effective media bulks almost coincides with the force between 
two-layered mirrors, VO2 on sapphire substrate. When the temperature overpasses the critical value, the entire bulk 
of the effective media becomes a metal, while in the case of two-layered mirror, the substrate remains an insulator. 
Consequently, at large distances the force between two effective media bulks is larger than the force between compound 
mirrors. The force obtained within the effective media calculation differs from the result of the two-layer calculation 
starting from the distance of several hundred nanometers (Fig.[6|). 

Finally in Fig. [7| we present the results for the VO2 mirror in front of the Gold mirror. The curves obey the same 
distance dependencies as in the case of equal mirrors. But the values of the reduction factor are increased thanks to 
the high reflectivity of gold. This setup is favorable from the experimental point of view. 

It is important to note that the dielectric function we use corresponds to thin film measurements. If the force is 
measured between true bulk VO2 mirrors, the results should be compared with the calculation that uses the bulk 
dielectric function [TT] ■ 



V. DISCUSSION AND CONCLUSION 



The obtained results can be understood in terms of the optical length of the layer material or the phase factor 
acquired by the field while propagating through the finite layer. 
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FIG. 5: Reduction factor of the Casimir force between two lOOnm (upper graph) or 4000nm slabs of Silicon with the carrier 
density modulated by laser irradiation in comparison with p-doped Silicon, n = 5 x 10^'' cm""^, tjp=0. 00184 eV, 7=0.00329 eV. 
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FIG. 6: Reduction factor of the Casimir force between two VO2 bulks (1), two 100 nm slabs of VO2 (2), and two 100 nm VO2 
layers on AI2O3 bulk substrate (3) below (a) and above (b) the critical temperature. 

As indicated in (U), the phase factor 6 is frequency-, wavevector- and thickness-dependent. Using the dielectric 
function of Silicon, which is typical for insulators, in expression (j4]) we find 



The phase factor ^ for metals or doped Si slabs modeled by a Drude or plasma model (7 = 0) evaluates to 



Let us now compare these phase factors at some constant fixed value for k, and D. Fig. [8] shows the phase acquired 
by the vacuum field while propagating through the intrinsic and doped silicon. The doping is described using either 
a Drude (dashed-double-dotted curve) or a plasma (dashed curve) model. 

The phase factors of intrinsic and weakly doped silicon are indistinguishable at low frequencies, provided that the 
doping is described by Drude model. Fig. [51 The phase factors acquired by the field in the silicon slab previously 
illuminated by laser, Fig.[9l coincide in Drude and plasma model descriptions down to the values of the ratio w/upe 
0.5. For lower frequencies the plasma model curve approaches the value \/2, and the Drude curve approaches zero. 
Though the difference between Drude and plasma curves reduces as the carrier density is increased, their behavior is 
qualitatively different in w, k — s- limit. It results in different predictions for the force at long distances. 

The largest contribution to the Casimir force comes from the frequencies around the characteristic frequency 
'-^ch ~ c/L. Large plate separation thus corresponds to small frequencies and small wavevectors k. At small frequencies 
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FIG. 7: Reduction factor of the Casimir force between VO2 and Au mirrors above and below the critical temperature. 
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FIG. 8: Phase factors acquired by the vacuum field when it propagates through the intrinsic and weakly p-doped silicon, 
cnjuji = 1, n = 5 • lO^cm-^, J^'^ = 1.84 • IQ-^eV, 7'"' = 3.29 • IQ-^eV, [l^. 
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FIG. 9: Phase factors acquired by the vacuum field when it propagates through the intrinsic and laser illuminated silicon, 
CKjm = 1, n = (2.0 ± 0.4) x 10^® cm"^ 



and small wave-vectors the bulk reflection coefRcients ^ tend to their static value 



lim II = ^^7== 12 



which is —0.55 for intrisic silicon, —0.52 for VO2 before the phase transition, and —1 for gold or any metal described 
by Drude or plasma model. Substituting these values into the integral ([2]), or in other words neglecting the dispersion 
at large separations of the mirrors, we get rough estimation for 77^' ~ 0.28 for two silicon bulks, r\^'^ ~ 0.25 for two 
VO2 bulks, and 77^"^* — 1 for two metallic ones. 

In [1^ we have shown that the effect of finite slab thickness manifests itself for Silicon at plate separations of the 
order L > c/luq, with ujq ~ 6.6 • 10^^ rad/s. 

The optical length of a silicon slab (jlOp at large plate separations tends to zero as Sins = D/cy^ e(iO) — 1 lu + 0{uj'^), 
while the bulk reflection coefhcient goes to its static value (IT2t . The numerator of ([4]) vanish. Consequently the 
reflection coeflicients of thin slabs at large plate separations vanish too. The force is considerably reduced. 

The optical length of a metallic slab (fTO|l at small frequencies and small transversal wave- vectors tends to a constant 
for the plasma model, S^^ — Dojp/c, while the bulk reflection coeflicients tend to —1 for both field polarizations. The 
reflection coefficient of a slab described by plasma model approaches its bulk value. The long distance limit is then 
not affected by the slab thickness except for very thin slabs. For gold with ujp = 9eV, two lOnm slabs separated by 
L = 10~^m yield up to the third decimal sign t] = F/Fc = 0.997. The reduction factor corresponding to 50 or lOOnm 
slabs is 77 = 0.999 which coincides with the one for two bulks at the same separation. 

For Drude model the optical length tends to zero as 5^^ = Duip/ {cy/^)^/uJ + 0(tJ^/^), while the bulk reflection 
coefficient for metals tends to p = — 1 + / u)p + 0{u}). The reflection coefficient is then r ^ 2p6/{l — p^{l — 2(5)). 
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Substituting here the low frequency expansions for p and 5 and confining ourselves to the lowest power of oj we get 

D 

The absolute value of the slab reflection coefficient (|13p is essentially smaller than one if -D < A. For gold with 
bjp = 9eT^, 7 = 0.035eV, the effective thickness is A = 1.7 x 10~^°m. Therefore in the long separation limit for slabs as 
thin as 20nm the reflection coefficients mount to —0.991! For two bulk mirrors separated by the distance L = 10~*m 
the reduction factor is 0.993. The slab reduction factor is lower than the bulk one at the same separation, but the 
effect is weak in comparison to dielectric mirrors. With the separation between plates L = 10~^m, the reduction 
factors are 0.922, 0.955, and 0.978 respectively for 10, 20 and 50 nm slabs. 

We conclude that it is enough to coat a silicon slab with a few nanometer layer of metal to suppress the effect of 
finite slab thickness. 

Though at low frequencies the phase factor acquired by the field in doped Si behaves similarly to the one in 
gold (Fig. [51 , the effective thickness A is much lower for gold than for doped Silicon where it varies from 7 x 10~'*m 
for the doping level of 10^^ cm~'^ to 64nm for ~ 10^° cm"'^. That is why we observe considerable reduction of the 
force already for 500 nm slabs. For 100 nm slab of doped silicon, n — 10^° cm~^, we get r ~ —0.61, and rj ^ 0.35. 
This rough estimation is in good agreement with our numerical result, rj ^ 0.38, see Fig. 01 

In section IV we presented the numerical results for VO2 which behaves as a semiconductor at temperatures below 
Tt = 340K and as a metal at higher temperatures. In the semiconductor state of this material the force between finite 
slabs reproduces the behavior characteristic for silicon. After the phase transition the slabs of VO2 are attracted as 
Drude metals. 

Of course all preceding results are based on a number of models, in particular the plasma and Drude models. As the 
issue of the precise description of finite conductivity for metals and semiconductors is still not satisfactorily settled, it 
would be very interesting to test the previously presented results experimentally. This might help to clarify the open 
problem of the temperature corrections to the Casimir force [3 EE, 113, 13 • 
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